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ABSTRACT 



Per certain assumptions, an analysis of multipla i^AD 
signals results in a reasonable estimate for the localiza- 
tion of a target relative to the SAD plaxforn. This is 
achieved by using selective approximations to linearize an 
initially nonlinear problem. The simulation ignores noise 
and requires an estimate of the magnitude of the target 
magnetic moment components. Results indicate that the best 
localization estimates are achieved when the platform is on 
cardinal headings, and when the target moment has a strong 
vertical component. 
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I. IMTRODOCTIDN 



A. BACKGROUND 

Ever since the operational introduction of one MAD 
(Magnetic Anomaly Detection) systaa in 1944, it has been a 
single sensor system limited to indicating detection only. 
The short ranges of the early versions allowed a general 
localization if a signal was detected. The use of multiple 
sensors should allow for a good estimate of the target’s 
location relative to the sensor platform, assuming that the 
arrangement of the sensors is such that the received signals 
are linearly independent of one another. Previous attempts 
to solve this problem have in general focused on the use of 
only two sensors, but did not develop a usable solution to 
operational completion. Jfork has been dene on this topic by 
Wynn et al. , [ Sef . 1], using a superconducting aradiometer 

array. The measured signals in their work were subjectsd to 
a novel signal processing technique which, in the labora- 
tory, was used to invert the dipole field equation to deter- 
mine both the position and the moment vectors of a dipole 
signal source. This thesis approaches the same proolem by 
simulating five standard total field measurements. As shown 
in Figure 1.1 , four magnetometers moving in parallel paths 
are used to obtain these measurements. The advantage of this 
approach compared to that of Wynn et al. is that it utilizes 
existing technology in conjunction with approximations, 
delineated in Chapter 3, in order to produce a usable 
result. 
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B. aPPBOACH 

The problem of localizing a target from any platform 
involves at least three unknowns: r^ , Cy , and , the compo- 
nents of the position vector from the platform to the 
target. In the case of using MAD, three additional unknowns 
niy, ffly, and aiy, the components of the target's magnetic 
moment, are also involved. The result is that there are six 
unknowns that must be solved for in order to localize a 
target using MAD. Since a normal sized magnetic targe* can 
be approximated by the magnetic field of a dipole [Hef. 2], 
it was decided to use the dipole equation as described in 
[Ref. 3] to generate the signals used in this simulation. 
Equation 1.1 describes the MAD signal, and will be developed 
from the dipole equation in chapter 2- Equation 1.1 is non- 

S (r) =- <m^+m=y+m^ /r3+3 (m^<- (re^^-yr^ej) /r s (eqn 1.1) 

linear in the unknowns. The position vector in equation 1. 1 
is defined as going from the target to the platform; this is 
a consequence of using the dipole equation, and will be 
explained later. 

The cartesian coordinates are defined in Figure 1.2 . 
The z axis is in the direction of platform motion, and the y 
axis is vertical down. The ccmpcnents of the earth's field 
unit vector, e^ , Sy , and e^ , are assumed to be known quanti- 
ties. Since it is very difficult to explicntly solve a set 
cf r.cn-linear equations in six unknowns, the problem will be 
simplified by making assumptions and dividing it into two 
special cases. The first case is a trivial on® in which the 
target's position is assumed known. For this case, equation 
1.1 is linear in the three unknowns a* , niy , and m^ . If the 
three measurements are linearly independent, then Cramer's 
Rule, [ Sef . U], may be used to solve for these unknowns. Th® 
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second case is the mere difficult one, and it is the prebiem 
investigated in this thesis. In this case, the assumption is 
made that the target's magnetic moment is known. That 
reduces equation 1.1 into a non-linear equation in only 
three unknowns, r^ , r^, and r^ . If rhis simplified version 
of equation 1.1 could be manipulated into a linear form by 
using several sensors and appropriate combinations of 
signals, the position of the target could easily be solved 
for. 



This thesis will develop an approximate solution to this 
problem by using combinations of the total field signals of 
an idealized multiple MAD system. 
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I 

Figure 1. 



Spatial Relatioaship of Sensor Positions. 
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I 
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Figure 1.2 Spatial Belationship of Platform and Target. 



II. THEORI 



The total field magnetometer measures the magnirude of 
the vector sum of the ambient magnetic field of the earth 
and rhs magnetic field of a magnetic dipole moment. The 
field of the dipole will used ro represent the magnetic 
field of the target. The magnitude of the vector sum is 
called the total field as given by 

K = + KlI 

represents the earth's magnetic field vector and 
represents the field due to a magnetic dipole. In the 
general case of interest here, the magnitude of the dipole 
field is much less than the magnitude of the earth's field. 
Under certain conditions, a useful approximation to the 
total field may be derived from the geometry shown in Figure 
2 . 1 . 




Figure 2.1 Vector Relationship of and B). . 

In this figure a very small vector is added to a very large 
one. The vectors of Figure 2.1 are not to scale. is to 

be considered several orders of magnitude smaller than . 
The application of basic geometric formulas for the sides of 
a triangle and the angle between sides yields 

3^2 = - 2b^B^cosJ^ 
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for the geooerry of Figure 2.1 . la terms cf e this equaiion 
becomes 

= BgZ + b^2 + 2b^BgCOse 

The 3^2 term can be factored out of the right hand side of 
the equation to produce 

= Bg2{i + (bi/Bg)2 f 2(bjL/Bg| C0S9} 

Taking rhe square root yields equation 2.1 . 

B^=Bg {[ 1 + (bj^/Bg) 2 + 2 (bj./B«)cos 0 ]V^} (eqn 2 . 1 ) 

Since by assumption bj^ is much less than B^, the rerm 
(bj_/Bg) 2 is much less than 2 (b^^/Bj ) cos9 and can be consid- 
ered negligible. This reduces equation 2.1 to 

Bt = S, {[1 + 2 (bjt/B*) cose] V^} (eqn 2 . 2 ) 

The general expansion formula given by 

(1+x)” = 1 + nx + n(n-1)x2/2! + Higher Order Terms 

can be applied to equation 2.1 . For the case here, n=1/2, 
and x= (bj^/Bg ) 2 2 (bj_/Bj) cosa . Substituting these values into 
the general expansion formula gives equation 2. 1 as 

B^ = Bg{1 + .5[{bji/Be)2 + 2 (bjL/Be ) cose ] 

- 1/8r(b^/3g)2 + (2bjL/3^) cose ]2 + H.O.T.} 

Squaring the appropriate term and multiplying through by the 
coefficients results in 

3^ = 3g (1 + .5[(b^/Bg)2] ♦ (b^/Be loose - l/Sf b^/Bg ] + 

- 4/8[ bj^ /Bg ]3cose - 4/8[ b^/3e ]2cos2 0 + H.O.T.) 

The cubed and higher order terms will be considered to be 
negligible. Dropping those terms and rearranging the 
remaining ones yields the form 
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= Be(1 ♦ (bj^/Be)cos0 - . 5 ) 200329 + .5(b^/Be)2} 

Usi.'ig the relation sinze + 00329 = 1, this raduoes rc 

B^ = Bg {1 + (bjL/Be)cos0 4- .5 (bjL /Be) 2 sin 2 8} 

After multiplying through by B^ , the seoond term car. be 
written in terms of the dot produot bj^_«Be,/3e . This is 
reduced to b^*® if fbe unit veotoc # is defined as Be/Ee» 
and the resulting equation is 

3^ = Be + b^*t + , 5 (bj_/Be ) 2sin20 

Sinoe bj_ is on the order of InT while 3^ is on the order of 
50,000 nT, the error term involving the sine is negligible. 
For the stated assumptions, equation 2. 1 oan be represented 
by 



B^ = Bg ♦ ^*e (eqn 2.3) 

In this thesis, 3^ will be assumed oonstann. In addition, 
this simulation will assume that no noise exists. Thus, the 
tonal field magnetometer will measure the anomaly in the 
earth's field defined by 

S (r) =5 

In this expression, bj^ represents the magnetic field of a 
dipole. bol is defined [Ref. 3] by the equation 

b^ (r) = -m/r2 + 3(m*r)r/r5 (aqn 2.4) 

m is the magnetic dipole moment in units of nTft^. r is the 
position vector to the field point and is measured in feet. 
If b^ is much less than 3^ , than aquation 2.3 and equation 
2.4 result in the equation 
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S (r) =- (m»e) /r 3 +3 (m «r) (r*9)/rs 



(egn 2.5) 



This equation is an alternative basis for the derivation of 
the Anderson's Functions, found in [Ref. 5]. The Anderson 
Functions are a commonlv used appro ximation for the signal 
received by a MAD sensor. The ccsaputar simulation in this 
thesis is based on equation 2.5 directly, since in is in a 
representation which is easier to use than the Anderson's 
function representation of a MAD signal. 



17 



III. NATORE OF THE PROBLEM 



This thssis will davelop and damonstrata an approxima- 
tion that gensratas a reasonable solarion of the targer's 
position relative ro a platform bearing several "sensors" to 
"measure" nhe magnetic signal as the platform makes a 
srraighr line encounter near the targe*. The assumpnions are 
as follows: 

1. ) The dipole equation is a valid representanicn of 
the magnenic field of the targe*. 

2. ) bj_ is several orders of magnitude smaller than . 

3. ) Measurements are done using only straighn line en- 
counters. 

4. ) The target's magnetic dipole momen* is a known 
vector . 

5. ) The earrh magnetic field vecror is known. 

6. ) The gradient of rhs earth's field is zero. 



7.) No magnetic ncise of 


any kind (environment. 


sensor , 


etc) exists. 






8.) The measurements 


are oriented along four 


lines 


which are oriented with 


respect to the flight 


path as 


shown in Figure 1.1 . 







Since the dipole equation is wrirten with the moment at the 
origin of a coordinate system, the calculations of the 
signal field values used in the simulation are made with the 
target located at the origin. As will be seen, the program 
performs a coordinate transformation in order to produce the 
output of target position relative to the platform. After 
this coordinate transformation, the origin is located at the 
sensor platform and moves with it. 

The five measurements are made along four separate 
flight paths, so only four different "sensors" are involved. 
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Sensor 4 could be replaced by inducing a rime delay in 
measurement of sensor 1. The signals generated at these 
measurement points are measured at each time step. The posi- 
tion vectors can all be described relative to the origin, 
located at the target, in terms of the position vector of 
sensor 1. The measured signals are represented by the 
following equations: 

S (f1)=- { (m,^+m^+m*) (r1) 2-3 ( nyc+nyy+n^z> ( x^+ y=y+z^ } /r 1 s 

S (r1) =- { + (r2) 2-3 ( m^ + ny; y-d]tn^) (x vC 1 i 

S (r5)=- { (m^+mye/m^ (r3) 2-3 ( mj[x-d ]+iy + n^) ([ x-d ]e^+ye + 2 ^ } /r3s 

S {r?i) =- { (m^e/ayy-m^) ( r4) z-3 ( nyc+ nyy+n^ z-d ]) (xe/y z-d } /r 4s 

S (r^) =- { (my-oiyy-m^y (r5) 2-3 ( m^x+d ]+a^y + e|Z) ([ x + d ]e^+ye+ze^ } /r 5 s 

For convenience, the following dummy variables are used in 
the computer program and are defined as: 



G = 


m^e^ + m^Sy ^-aije^ 


VI 


- 


3 m ^ “ G 


V2 




3 (m^ey + raye^ ) 


73 


= 


3mySy - G 


7 4 


= 


3 (fflye^ + m,ey ) 


75 


= 


3m,e^ - G 


76 


= 


3(m,e^ - m^e^) 


S 1 




S (f1) 


S2 


= 


S (f^) 


S3 


= 


S(r^) 


S4 


rr 


s (ffi) 


35 




S(r^) 



The use of these variables reduces the above equations to 
the following form, which are the defining equations for the 
signals in the simulation. 

S 1= (V1x2 + V2xy+ V3y2 + 74yz<-V5 z2+ 76xz» /r is 
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S2= (Vlx2 + V2x7-V2xd + V3y2-2y 3yd+V3d2 + 74yz-V4dz+V5z2 + V6xz) /r 2 s 
S3= (yix2-2V1xd+71d2 +V2xy-7 2dy+73y2<- V4yz+V5z2+V6xz-76dz) /r35 
S4= (71x2+V2xy+73y2 + 74yz-74yd+75z2-2V5zd+75d2 +76xz-76xd) /r4s 
S5= (71x2+271xd + 71d2 +V2xy+7 2dy+73y2 + 7 4yz + 75z2 +V6xz + 76dz) /rSs 

The unknowns of interest are x, y, and z, which are related 



to the position vectors as follows 


: 






Fl = X"! + yg* + z1c 




(eqn 


3.1) 


r^ = xt + (y - d ) g* + z1c 




(eqn 


3.2) 


ri = (x-d)t + yT + 2^ 




(aqn 


3.3) 


= x“? + yg* + (z-d)1c 




(eqn 


3.4) 


ro = (x+d)“^ + yg“ + z1c 




(eqn 


3.5) 


Tc linearize the problem it 


will 


be assumed tha 


^ the 


denominators of the equations for 
given above are all equal to r, i. 
r = r1 = r2 = r3 •= r4 = r5 


SI , 

3. f 


S2, 33, 34, a 


nd S5 


The equations for SI, S2, S3, 34, 


and 


S5 with the denomina- 


tors equal can then be subtracted 
the following: 


from 


one another to 


yield 



S2-S3={(72-27 l)x+ (273-72) y+ (74-76) z+ (71-73) d} (-d/rS) 

52- S5={(72+271) x+ (273+72 )y + (74+75) z+ (71-73)1} (-d/rS) 

53- S1= {27lx + 72y + 78z-7 Id} (-d/rS) 

S5-S 1= {- 271 x-72y-76z-71d) (-d/rS) 

S3-Sa={(271-76) x+ (72-74) y+ (76-275) z+ (75-71 )d} (-d/rS) 
S5-S4= {- (271 + 76) x-(V2 + 74) y- (76+275) z+ (75-71) d} (-d/rS) 

Sines the denominators are assumed to be identical, if these 
six equations are divided in pairs the coefficients 
involving the non- linear r^ term and the quantity d are 
canceled out. This results in the following equations that 



define the new variables A1 , A2, and A3 as 






A1 = (S2-S3) / (S2-S5) 


(eqn 


3.5) 


A2 = (S3-S1) / (S5-S 1) 


(eqn 


3 .7) 


A3 = (S3-S4) / (S5-S4) 


(eqn 


3.3) 
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Defining the follnwing variablas for convenianca. 

Cl = A1 (V2+2V1) - (V2-2V1 ) 

C2 = Al(2V3+72) - (2V3-72) 

C3 = A1 (V4+76) - (V4-76) 

K1 = [ (71-73) - A 1 (71-73) ]d 

C4 = A2 (271) + 27 1 

C5 = A2(72) + 72 

C6 = A2 (76) + 76 

K2 = C VI - A2 (71) ]d 

C7 = A3 (76+271) - (76-271 ) 

C8 = A3 (74+72) - (74-72) 

C9 = A3 (275+76) - (275-76) 

K3 = [ (71-75) - A3 (71-75) ]d 



equations 3,6 - 3,3 can be written in rhe folloving form as 



C 1 X 


+ C2y 


+ C3z = 


K1 


(eqn 3 


.9) 


C4x 


+ C5y 


+ C6z = 


K2 


(eqn 3. 


10) 


C7x 


+ C8y 


+ C9z = 


K3 


(eqn 3. 


11) 


Thes e 


t hr ee 


equa tion 


s are linear in terms cf x. 


y, and z. 


In 



addition, they are linearly indepan ienn as shown by using 
the program in Appendix B, A sample ounpun demonstrating the 
linear independence is shown in Appendix C, Since nhey are 
linearly independent, Cramer's Rule can be used to solve 
equations 3,9 - 3,11 for x, y, and z. The computer program 
in Appendix A was written in terms of these equations. The 
program assumes that the magnitude of the target's magnetic 
moment is 5x19® nTft^, A representative moment was selected 
by using the values listed by Fromm in [Sef. 2], who lists a 
moment of 10® to 2x10® cgs units for a submarine. One cgs 
unit is approximately 3.35 nlft^. The value used in the 
program is a rounded out average of Fromm's values. 
Computations are done using the standard convention with the 
target dipole at the origin. However, as indicated earlier. 
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the output is corrected by a coordinate transfer mati on so 
that the platform is at the origin, and the target position 
is given with respect to the platform. 
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17. DISCOS SID H 



A. GENERAL 

The program in Appendix B was written to check whether 
the results generated by the approximation equations were 
reasonably close to the actual coor dinar es of the rarget. 
The program was modified inro the fora shown in Appendix A 
to graphically provide a comparison between the calculated 
and the actual position coordinates. Initial parameters for 
the knowns in the equations were a platform heading of 30° 
magnetic, an earth field vector of 70<3 down from the hori- 
zontal, an d a target dipole moment vector with a vertical 



component 


500 down from the 


horizontal 


an d 


a horizontal 


component 


oriented 355o from 


magnetic nor 


th. 


These initial 


parameters 


have no special si 


gnifirance , 


and 





sent a convenient starting point. As discussed earlier, the 
magnitude of the target dipole moment was set at the value 
of 5X10® nTft®. The spacing between measurement positions, 
d, was taken as 50 feet, and the platform passed directly 
overhead the target at 1000 feet. This and all subsequent 
runs were started at 5000 feet prior to CPA (closes* point 



of approach) with 


the platform aiva 


nc ing 50 


feet per 


• ime 


step. One position 


calculation was 


CO mpletad 


at each 


time 


step, the simulation 


ceasing when the 


plat for m 


was 5000 


feet 



past CPA. For interpretation purposes, calculated target 
positions are shown as the sensor platform moves 10,000 feet 
past the target in a straight line encounter, with CPA 
occurring at the 5000 foot point. Neither the CPA nor the 
target’s location are in any way considered known in the 
localization process. However, as the area of interest is 
the behavior of the simulation close to CPA, the closest 
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AS ACTUAL 



Figure 4.1 Initial Run, x position, 

point of approach is always centarsd in tha fcllcwing 
figures. 



24 



g =aa|l ^B jLj jj ,^ipi|p '• 

OiP. ' jia 






lllllll■■lllJI>lrl^l» ltill ■ %>>»..*i-'-»:»te|llliillllfcg;fiwl^^ ,htfc uVr 






MWIi,.-^^ 









'■iNfjw:!:: 



bk Ti niiiiiiitiii HI jpoiBHHni 
M’ i-,„— T 

m IBl :v'r , , 1.*, 






mi m 



«» W7I 



wmm 








YPOSmON 




OSTANCE ALONG TRACK > 

LECCNO 

o s CALCULATED 
asACTUAL 



Figure 4.2 Initial Run, y position. 



From the orientation 



of the coordinate system shown in 



Figure 1.2 , z is in the direction of flight, y is straight 
down, and x is off the right side of the platform. Since 
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Figure 4.3 Initial Hun, z position. 

this thesis investigates straight line encounters only, 
x(actual) and y (actual) coordinates will remain const 
throughout any run for the coordinate orientation used. 



the 

ant 

The 
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CORRECTED Y POSmOH 




LICCNO 

o s CCRR. CALCULATED 
asACTUAL 



Figure 4.4 Initial Run, corrected y position. 

five maasursmen t s will tharefcrs be made wi*.a their respec- 
tive X and y coordinates as constants. The z coordinate will 
vary. This means that for the initial pass described. 
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x(actual) would be constant at zero, y (actual) would be 
constant at 1300 feet, while z (actual) would progress from 
5000 feet to -5000 feet. Figures 4. 1 - 4.3 show a compar- 
ison of the actual and calculated target positions as -he 
platform completes the initial run described earlier. It 
will be noted that the calculated and actual x, y, and z 
positions are almost identical at the start of the run. This 
is a coincidence. "Crossover” points such as this occur for 
all three coordinates simultaneously during a run, falling 
closer or further away from CPA, depending upon the choice 
of parameters. 

In this initial run, it is obvious that the "worst" 
estimates occur for the y coordinate. This is less signifi- 
cant than it seems at first glancs, however, as the y coor- 
dinate is the vertical separation between the platform and 
the target. If the target is submerged, this value is 
unknown, but the altitude of the platform is known. The 
target can be no closer to the platform than thrs altitude, 
since the target is physically unable to rise above the 
ocean surface. As an example, if a target dapth of 200 f set 
is assumed, the initial y separation of 1000 feet forces the 
altitude in this simulation to be 830 feet. Using this value 
as a bound on y , it can be seen in Figure 4.4 that the y 
calculated position is somewhat improved. Accordingly, it 
seems that the particular arrangements of sensors shown in 
Figure 1.1 does show promise in the problem of utilizing 
this geometry of Multiple MAD for accurate localization 
purposes. Accordingly, variations of parameters were tried 
in order to determine what effect ohanging them would have. 
In all cases, parameters not specifically mentioned in a 
section will retain the values as for the initial run 
described above. The singularity-like behavior occaring at 
various places on the graphs are a result of the determinant 
of the denominator in the Cramer's Rule subroutine having a 
value close to zero. 



28 



B. HEADING CHANGES 



In -chis saction, the localization process will be inves- 
tigated as a function of heading. Compared to rhe initial 
conditions described above, the crossover points were found 
to be moved further cur past 5003 faat prior to CPA on a 
heading of Oo. While a specific example of this is not 
shown, the behavior of the crossover poinrs for the plarform 
headings of 30<> and 60° may be observed by comparing Figures 
4.1 - 4.4 to the respective portions of Figure 4.5 . The x 
position estimates were better for a North-South heading 
than for any other platform heading, although the estimares 
on East-West headings were also very good. For a 90° heading 
there are two crossover points located equidistant at 1200 
feet from CPA. Calculated positions on all headings are 
generally good, excepting those around 60®, where the x 
estimates shows a large discrepancy at distances greater 
than 3000 feet prior to CPA. The y estimates are not signif- 
icantly affected, but the z estimates are good only within ± 
1200 feet of CPA. Nevertheless, even in this worst-case 
situation, reasonable calculated coordinates are available 
for localization along a significant portion of the track. 

As shown in Figures 4.5 - 4.8, headings in the quadrant 
Qo - 900 reasonably describe what happens on any heading, as 
the graphs for analogous headinas (50o, 120®, 24Qo, and 30Qo 
in this example) have similar shapes, although the orienta- 
tion may differ between quadrants. These similarities are 
apparently dua to the fact that analogous heading in quad- 
rants 1 and 3 and quadrants 2 and 4 are in fact reciprocal 
headings, and so one would expect ohs graphs to maintain the 
same shape; however, the ending point of one would represent 
the starting point of the other. The relationship between 
adjacent quadrants probably exists for similar reasons. 



29 



C. TARGET DIPOLE HOtJENT CHANGES 



In this section, the horizontal component of m was 
oriented in the 0° (magnetic North) direction vice 355o, and 
the vertical component of m was varied. The best positioning 
was obtained for m oriented vertically up,i.e., no hori- 
zontal component, and an example for that case is shown in 
figure 4.9 . Poor positioning was obtained for the vertical 
component of a pointing between the horizontal and 45° up. 
The worst case shown is Figure 4.10, where m is oriented 
horizontally pointing due north, w'ith the moment oriented at 
an angle below the horizontal, reasonable positioning esti- 
mates were obtained. There was less variability in posi- 
tioning for moments pointing below the horizontal than for 
those pointing above the horizontal. 

The results for variations in the orientation of the 
target dipole moment suggests than the positioning estimates 
are less sensitive to moment orientation changes in the 
horizontal plana as compared to tae verrical plane. For 
example, by keeping the vertical componenr of m constant at 
an angle of down while varying the horizontal component 

direction it was found than the horizontal component 
affected the positioning much less than variations in the 
vertical. Illustrated in Figure 4.11, the x estimates of 
position improve as the horizontal coraponenr approaches the 
plarfcrm heading, while the y esrimares improve as -^he hori- 
zontal component approaches alignment with magnetic North. 
The latter case is not shown. 

D. CHANGES IN SENSOR SPACING 

Changing the inter-sensor spacing had at mosr a weak 
effect upon the position outpur, even wirh spacing as short 
as 5 feet or as long as 200 feet. This would be expected in 
a computer simulation where sufficient precision is 
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E. ALTITODE CHANGES 

For a straight overhead pass, positioning inforxat ion 
appears best if a pass with a 500 to 1000 foot separation 
between the target and the platforx is made at CPA. At 
greater altitudes, the accuracy suffers slightly with the 
increase in altitude, while at lesser ranges, usable posi- 
tioning information is obtained only relatively close to 
CPA. For example, for a vertical separation of 400 feet, 
positioning data was very inaccurate at slant ranges greater 
than 1000 feet. 

F. NON-OVEBHEAD PASSES 

This program gives poor results- for passes that are nc^ 
directly over or directly off to one side of the target. The 
worst situation observed in any variation is illustrated in 
Figure 4.12 . To generate this figure, a pass was made such 
that the actual x coordinate was approximately equal to the 
actual y coordinate. The best positioning was achieved by 
overhead passes, but low off to the side passes can yield 
good results as well. 



G. AVERAGES 

The averages resulting from a number of runs are 
presented in Table I . These selected averages, covering 
all of the variations discussed in this chapter, show little 
consistency. This is probably due in part to the existence 
of the singularities. The variability is also a function of 
the limits of travel for the runs. Nevertheless, for most 
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Figure 4.10 la Horizontal and Oriented to Hagnetic North. 
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Figurs 4.12 Target at 45<> depres. angle rel. to platfora at CPA 
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cases, they do show reasonable values, and by thamseives io 
provide limited localization data. The problem for inflight 
analysis would be to determine the bounds of a run in rhe 
absence of any other target positioning data. 

H. IMPOETANCB OF C0NSIST3NCI 

An important observation resulting from the analysis of 
the graphs, including those not shown as examples in this 
chapter, is that the rate of change of the calculated coor- 
dinates gives some indication of the validity of the posi- 
tioning data. If the x and y grapas have a slope close to 
zero, and if the z graph has a slope of -45<J for the calcu- 
lated coordinate, then positioning data is likely to be 
comparatively accurate. If all three graphs fulfill this 
condition simultaneously, then positioning data appears to 
be usable. 

I. IMPLEMENTATION 

The possible implementation of the process described in 
tnis thesis would be to assume a magnetic moment for the 
target, then run the signal data through the processing to 
generate positioning data, and then compare these results to 
those obtained by using the assumed moment to generate a 
signal via this simulation. If the agreement between results 
were good, cne cculd assume that one had assumed the correct 
target moment; if not, one would assume a different moment, 
and go through the process again, until one had a reasonable 
target moment, which would be valid at least until the 
target changed course. Osing a computer these comparisons 
cculd be made reasonably quickly, so this process might have 
operational possibilities. 



V* COHCLOSIOH 



For the assumpticns made, the positioning estimates are 
reasonable. Also, the aysrages taken of the three coordi- 
nates as calculated throughout the various runs show that 
the X calculation averages are in general at worst within 
500 feet of the actual x value, while the worst y calcula- 
tion averages tend to be within 1000 feet of the actual y 
value. The best calculated x and y values were within 100 
and 200 feet, respectively, of the true values. 

Trends observed were as follows: 

1. It is best to fly the sensor platform directly over 
the target, or alternatively to have a 'distant' 
(between 1000 and 2000 foot) CPA off to one side while 
minimizing the target to platform separation in the y 
direction. The worst situation is when the target is at 
a point that, at CPA, has a depression angle of 45° with 
respect to the platform. 

2. Best vertical separations are those that will yield a 
500 tc 1000 foot vertical distance between target and 
platform. Short vertical separations should be avoided 
if good localization is desired. 

3. Best results are obtained if headings are northerly 

or southerly. Poor results are achieved wi'-h headings 
close to 60°, 1200, 240°, and 300°. Sast-tfest headings, 

while not as good as North-South, do provide good posi- 
tioning data close to CPA. 

4. Good results may be expected if there is a large 
vertical moment, especially if it is vertically up. 
Poor positioning results are obtained with this simula- 
tion if the moment has a large horizontal component. 
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5. Results are slightly better if the rarget momenr is 
lined up with either the axis of flight or with the 
earth's field, although this effect is ofren doniinarsd 
by the more gross effects caused by rhe other 'variables. 

A suggested approach for further study is to inject 
noise into the simulation in order to see how robust it is 
in the accuracy of the position calculations. If the 
results remain reasonable, this approximation or an altered 
form of it may one day have real-world applicability. 
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APPENDIX A 



COMPUTER PROGRAM FOR PLOT OUTPUT 
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APPENDIX B 



COMPUTER PROGRAM FOR NUMERICAL OUTPUT 
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APPENDIX C 

EXAMPLE OOTPDT OF APPENDIX B 



An axampls 
shown below. 



of the output of tha 
The total output fo 



program in 
a single 



Appendix 3 
time steo 



is 

is 



presented. 



SENSOR 

number 

1 

2 

3 

A 

5 



EQNS 
USED 
1 2 3 



SLhNT 

RANGE 

5C50.0 

506C.1 

5050.2 

509S.0 

5050.2 

ECLATION 

S2-S3 

52- S5 

53- S1 
S5-S1 
S3-SA 
S5-SA 



ACTUAL 





SENS ED 
SIGNAL 
0.60540-04 
0.1 5960-C3 
0.2682 D-CA 
0.5 9400-04 
0.93660-CA 



CALCULATEO 
SIGNAL 
0.6054E-04 
0.1 612E-03 
0.2683E-04 
0.4135E-04 
0.9 3&8E-04 



TRUE 

0.1327C-03 
0.6590C-04 
-. 3372C-04 
0. 5212C-04 
-. 1258C-04 
0. 5A26C-04 



APPROXIMATE 
0.1343E-03 
0.o7e8E-04 
-. 33 7 1 £-04 
0.3314E-04 
- .1452E-04 
0 .523 3E-04 



X Y 

POS POS 

0.0 - 1000.00 



z 

PCS 

-4950. 00 



The first fi/a rows of data above are for the appropriate 
‘'sensor" as indicated, the slant range from that sensor to 
the target, the received signal as generated by equation 2.5 
and the signal that would be generated using the approxima- 
tion described in Chapter 3 to linearize the problem. The 
next six rows show the results of taking the differences 
between the received and calculated (for comparison purposes 
only) signals, listed under the headings of TRUE and 
APPROXIMATE, respectively. The final two lines shew the 
actual coordinates of the platform relative to the target 
(the coordinate transformation has not been done at this 
point) as compared to the calculated position. These are the 
points used by the program in Appendix A to plot the graphs 
after the coordinate transformation is completed. 
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The seconi secricn of the output is showr. below 



EiJATIQIJ 

NJMBER 



CONSTANT 



X 

COEFFICIcNT 
-0. 30D+02 
0. lAJ+UC 
-0 • 21D+01 



COEFFICIENT CCEFFiCicNT 
-0.920+01 0.200+01 

0. iSO+OO u. 540-0x 

—0.190+02 0.5oU+0i 



y 



z 



1 

2 

3 



0. 5 00+ 03 
-0. 33o+03 
-0. 6 o 0+ 02 




-0.i3£+01 
— 0. 5 8 E+ 0 1 
0.44E+01 



'0.22t + 03 -0. 71E+02 0.3oc + 02 

0.14c+02 0.t9E+00 0.5oE+00 

0.63c— O a -0. 69c— 02 0.15c— 01 



The first rhree lines of the ouxput represent the coeffi- 
cients of the equations 3,6 - 3.8 in the form of 

CONSTANT = X COEPF(X) + Y COEFF(y) *■ Z COEFF(z). In order to 
determine the linear independence of these equations, the 
ratios of the coefficients are taJcen in pairs. Since it may 
be seen from the lasx three lines of rhe outpux that the 
ratios of the constants and coefficients differ signifi- 
cantly for all three equations, it ran be concluded that the 
equations are linearly independent. This holds true 
throughout any given run, which may be verified by printing 
out all time steps for that run. 
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AP PENDIX D 
CRAMER'S RULE 



A. GENERAL 



Cramer's Buis, a method for solving systems 
equations in tvo or more unknowns, states that sa 
can be expressed as the ratio of two datermina 
general form of a solution is shown below. 



Ax+By+Cz=o 
Bx+Fy+Gz=H 
I x+ Jy+Kz=L 



D 



ABC 
£ F G 
I J K 



X 



0 B C 
H F G 
L J K 



/ D 



y 



A 0 C 
Z H G 
ILK 



/ D 



z 



ABO 
B F H 
I J L 



/ D 



B. EXAMPLE 



The following sxampla was taken from [Hef, 4], 
3x-2y+2z=7 
x+ y+ z=6 
2x- y-2z=2 



of linear 
h unknown 
es. The 



Sol ve : 
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Ler D be the denominator. 





1 


3 


-2 


2 




D = 


( 


1 


1 


1 


= -17 




1 


2 


- 1 


-2 






1 


7 


-2 


2 1 


1 




1 


6 


1 


1 1 


1 -51 


X = 


2 


2 


-1 


-2 


1 = :i7 








D 








1 


3 


7 


2 






1 


1 


6 


1 


-34 


Y - 


1 


2 


2 


-2 


= 












-17 








D 








1 


3 


-2 


7 1 


1 




1 


1 


1 


6 ! 


1 _ -17 


z = 




2 


-1 


2 i 














- ~ -17 



D 



CHECK: 

3x-2y+2z=7 
3 (3)-2 (2) +2(1) =7 
9-4+2=7 
7 = 7 

x+y+z=6 
3+2+1+6 
6 = 6 

2x-y-2z=2 

2(3)-2-2(1)=2 

5-2-2=2 

2=2 



C. SOBEOOTINB "CR1MEH*S ROLE” CHECK 

Running the values in Section 3 through the subrourinas 
of the programs in Appendix A and Appendix 3 yielded nhe 
same results as the example. 



58 



LIST OF RBFEEBSCES 



Hynn.W.a., Frahm ,C. P. , 
Weilhoner, J. , Mynn,M.J, 
qradioms t e r / a a qnq roae ter 
proce ssi ng •e^nlgue, 
Hagne'ics, 



Carroil,P. J. , Clark, R.H., 
Advanced superc c r.d a ct iv e 
arrays an3 a novsT^ signal 
-p77UT-7U77 ■ Tfans7 



Froinm,W. S, in Advances in Electronics, V. 4 , edited by 
Marten, L., p . 2507 '5ca‘3emTc ^ress,”T752 



Reitz,J. R., Milford, F.J. , 2- 

Elect rom a gne t ic Theory, ?.164, lidaison-Wesley 

FiT5IIs^ng!:o77 Inc77”l9 67 



Shapiro, M . S. , Mathe matic c Encyc lope dia. p. 49, 
Doubleday 5 Co moan y7 “Inc. , T977 ~ “ 



Anderson , J . E. , NADC Report N ADC-EL-47-50 , Ma o ne ti c 

Ai rbo rne Det e ction Frequency Fesoonses, 1949 ~ “ 



BIBLIOGBAPHY 



BOEING Aerosoace Company Proprietary Data, N A 7Y/ BOSING 
Development of a D ual a AD Loc a liza tron Sy ste m, v em'Ee? 



BOEING Aerospace Company, Dual MAD Processinc Stud^, 1 

A dvanc ed MEG Pro gra m,~ 1^ov = m59r~lBTg " “ ~ 

BOEING Aerospace Company, Dual MAD Processi^ Studv, Tas k 2 

Report , Techni q ues A nalysj gT T UeSemSe? IVlt ~ 



National Defense Research Commitnae Division 6 Technical 
Report, Magnet ic Ai rborn e D etecrion P rogra m , V.5, 1946 



Naval Air Development Center Technical Memo 
NADC-30 1 2 : ARO, Extrema, symmetry, and zeroes of HAD s 
and t_^ir a cpl i caf Ion to a^IPl oca liza^ion , ~7U Dctooer 



r an dum 
ignals 
TBB2 — 



TEXAS INSTRUMENTS INCORPORATED, Analysis o^ Transverse 

Gra diomete r MAD usi^ Dual Sca lar gag ner oms ters ?aff~~TT 
ITefecrro n ~goc[e I ancT~ 7err ormance, ITcGf egor , D. , July 1976’ “ 

TEXAS INSTRUMENTS INCORPORATED, Analysis of Transverse 

Gradicmerer MAD using Dual Scalar Magneno me t efs, Part ' TIT 
Tac?rcax ~ Param e tef '~ '5xTra crr 9P , !TcGregof , D. , January T978 

TEXAS INSTRUMENTS INCORPORATED, Flight Tesr Plan Dual M^, 

McGregor,D., undated " " “ “ 



60 



INITIAL DISTBIBOTIOH LIST 



Nc. Copiss 

1. Fleet Tactical Library 2 

P.O.Box 104 2 

Silver Springs, Md. 20910 

2. Defense rachnical Inforniation Center 2 

Cameron Station 

Alexandria, Virginia 22314 

3. Library, Code 0 142 2 

Naval Postgraduate School 

Monterey, Ca. 9 3943 

4. Dr. a. N. Forrest, code 7 1 2 

Chairman ASW Academic Group 

Naval Postgraduate School 
Monterey, Ca. 93943 

5. Dr. Otto Heinz, Code 6 1Hz 2 

Department of Physics 

Naval Postgraduate School 
Monterey, Ca. 93943 

6. Dr. Andrew R. Ochadlick, Jr., Code 610c 2 

Department of Physics 

Naval Postgraduate School 
Monterey, Ca. 93943 

7. LCDR C. Burmaster, USN, Code 61Zr 1 

Department of Physics 

Naval Postgraduate School 
Monterey, Ca. 93943 

8. LT Wolf H. Bocic, OSN 2 

ASWOC Sigonella 

FPO New ?or)c, N.Y. 095 23 

9. Mr. John Shannon 1 

Naval Air Develooment center (Code 30 12) 

Warminster, Pa. 18974 

10. Mr. Edward Yannuzzi 1 

Naval Air Development Canter (Code 30) 

Warminster, Pa. 18974 

11. Chief of Naval Research 
Deoartment of the Navy 

Attn: John G. Heacocic (Code 425 GG) 1 

800 North Quincv Street 
Arlington, Va. 22217 

12. Mr. William Andahazy 1 

Naval Ship Research & Development Center 
Annatolis Laboratory 

Annapolis, Md. 21402 



6 1 



2025G0 



j Thesis 



Thesis 

B58895 
c .1 



Bock 

A linear approxima- 
tion of the source _ 
position using multi- 
ple MAD. 




Thesis 



SFP SO 



3 7 0 5 3 



20 



o 




B58895 Bock 

^^2. A linear approxima- 

tion of the source 
position using multi- 
ple MAD. 



cn 



